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Abstract. - The electronic structure of p-type doped Bi2Te3 is studied by angle resolved photoe- 
mission spectroscopy (ARPES) to experimentally confirm the mechanism responsible for the high 
thermoelectric figure of merit. Our ARPES study shows that the band edges are located off the T-Z 
line in the Brillouin zone, which provides direct observation that the spin-orbit interaction is a key 
factor to understand the electronic structure and the corresponding thermoelectric properties of 
Bi2Te3. Successive time dependent ARPES measurement also reveals that the electron- like bands 
crossing Ef near the F point are formed in an hour after cleaving the crystals. We interpret these 
as surface states induced by surface band bending, possibly due to quintuple inter-layer distance 
change of Bi2Te3. 



It is of great interest to find a good thermoelectric (TE) material competitive to conven- 
tional compressor-based active cooling/power generation systems. Desirable features such as 
lack of noise or emission of polluting gases, space saving, high reliability and semi-permanent 
life make them very attractive though to date their main shortcoming has been their low 
efficiency [1]. Appropriately doped Bi2Te3 has been one of the most competitive candidate 
materials for fulfilling the above purpose since its high efficiency was discovered 70 years 
ago [2]. Since that time, much effort has been expended to improve its TE properties by 
doping and alloying, but little progress was made until the reports for P^-^Shr Te3- y Se y al- 
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loys [3]. Even though higher efficiency at room temperature was reported in Bi2Te3/Sb2Te3 
superlattice [4], it is very important to completely understand the electronic structure of 
Bi 2 Te 3 in order to establish a reliable guideline for finding or synthesizing better TE com- 
pounds. 

Narrow gap semiconductors composed of covalently bonded heavy elements have gen- 
erally been considered as prominent TE materials because the gap and covalent bonding 
(heavy elements) have a tendency to increase (decrease) the electron entropy and mobil- 
ity (the thermal conductivity), respectively [5]. In addition to these general tendencies, a 
few system-dependent features have also been pointed out for Bi2Te3; highly anisotropic 
(pseudo) two-dimensional structure [6]; the six- valley electron (hole) pockets in n(p)-type 
doped Bi 2 Te3 [7-10]. Among these features, the six- valley electron and hole pockets in the 
electronic structure of Bi2Te3 drew considerable attention. Since the multivallcy conduction 
band structure of Bi2Te3 was inferred from magnetoresistance and Faraday rotation experi- 
ments [7,8], the positions and angles of the hole and electron pockets have been investigated 
by de Haas-van Alphen and Shubnikov-dc Haas experiments [9,10]. However, a controversy 
about the six- valley model arose in theoretical studies. The multi- valley structure had not 
been perfectly reproduced by the band calculations at an early stage [11-13]. This con- 
troversy was not settled until spin-orbit interactions were included in the band calculation 
scheme. Even in this case, not all the methods including spin-orbit interaction give the same 
results. Early calculations gave six valence band maxima (VBM) and two conduction band 
minima (CBM) [14-16], but recent ab initio calculation studies give six extrema for both 
the valence bands and the conduction bands, showing that the spin-orbit interaction is a key 
factor for understanding the electronic structure and the TE properties of Bi 2 Te 3 [17-20]. 
However, the issue is not clearly resolved yet. Although the six-valley model has been pro- 
posed for forty years and many theoretical works have been undertaken to explain it, the 
model itself is based on only indirect experimental facts; no direct evidence has not been 
provided so far. Establishing direct evidence is still crucial for confirming whether the spin- 
orbit interaction is really a dominant factor in determining the electronic structure of Bi 6p 
orbitals hybridized with Te 5p orbitals. 

Angle-resolved photocmission spectroscopy (ARPES) is the most direct method for inves- 
tigating the electronic band structure of solids. Photoemission process in a solid corresponds 
to single-particle excitation from an occupied state to unoccupied one, so the photoelectron 
contains information of the energy and the momentum inside the solid, i.e. the electronic 
structure of the solid. Owing to this direct connection between the observables and the 
objectives, ARPES plays an important role in researching hot-debated materials. To date, 
we are aware of only one ARPES experiment on Bi 2 Te3 that has been reported in the lit- 
erature [21]. The angular resolution of the spectra from that experiment was insufficient to 
address the issue cited above. In this paper, we present a direct observation of the six-valley 
model in p-type doped Bi 2 Te3 by ARPES. Our ARPES measurements with a high angular 
resolution reveal that the valence band maxima (VBM) are located off the T-Z line in the 
Brillouin zone (BZ) as predicted in the ab initio band calculation including spin-orbit in- 
teraction [17-20]. Further, we find that surface states are formed near T point in an hour 
after cleaving the sample in situ, which has obstructed direct observation of the correct bulk 
electronic structure of Bi 2 Te3. 

Two kinds of Bi 2 Tc 3 single crystals, Bi 1% excess (Bi-1%) and Te 4% excess (Tc-4%), 
were grown by melting Bi and Te in a mixture with the molar ratio and slowly cooling. 
Resistivity measurements using the standard four probe method show both of samples (Bi- 
1% and Te-4%) are metallic. Sccbeck coefficient and Hall coefficient measurements show 
the samples are p-type. For ARPES experiments, the samples were attached to a cryostat 
and introduced into our analysis chamber in ultra-high vacuum environment of 5xl0~ n 
Torr. A shiny cleaved surface was obtained by the top post method in situ, and the sample 
orientation was checked by low energy electron diffraction (LEED). Low photon energy 
measurements (hu> ~20 eV) were performed on undulator beam line U13UB at the National 
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Fig. 1: (Color on-line) (a) Crystal structure of Bi2Te3. The red rhomb ohedron represents a rhombo- 
hedral unit cell, (b) Rhombohedral Brillouin zone of Bi2Te3 . The scale of the drawing is exaggerated 
along rZ for the purpose of clarity. The six orange spheres represent the hole pockets in the six val- 
ley model. The small blue segment containing one hole pocket represents the measured k position 
of the valence band ARPES spectra. 

Synchrotron Light Source with a Scienta SES-2002 electron spectrometer. High photon 
energy measurements (Tiuj ~121 eV) were done on undulator beam line U10 at the Pohang 
Light Source. The instrumental energy resolution and momentum resolution were set to 
20 meV (30 meV) and ±0.1° for the low (high) photon energy, respectively. The reported 
spectra here, if not otherwise specified, were obtained at 20 K within 15 minutes of cleaving 
the samples. 

The crystal structure of Bi2Te3 is rhombohedral with the space group D^ d (Rim). The 
corresponding Brilloune zone (BZ) with high symmetry points is depicted in Fig. [TJ Along 
the rhombohedral [111] direction, so-called quintuple layers (Te(l)-Bi-Te(2)-Bi-Te(l)) are 
stacked, forming intralayer Bi 6p-Te bp covalent bonding and interlayer van der Waals bond- 
ing. The measured k positions in the BZ are depicted as a blue segment in Fig. [IJb). In 
principle, the k z value can be determined if we use photon energy dependence of the pho- 
toemission spectra. Unfortunately, due to the surface effects, discussed below, the band 
dispersion along the k z direction is very weak. Because of this attempts at determining 
the position of k z by changing the photon energy were not successful. Instead, we use the 
constant inner potential approximation with the inner potential, Vq= 10 eV, to determine 
the value of k z [21,22]. The normal emission spectra taken with 7iu;=24.0 eV then approxi- 
mately corresponds to the Z point in the BZ, and so we chose ?Lw=23.0 eV to search for the 
valence band maxima in this study. 

ARPES images of the valence bands of Te-4% are shown in Fig. HJa). The valence 
band structure of Bi-1% (not shown) was measured and found to be very similar to that 
of Te-4%. For comparison, the bands calculated by the local density approximation (LDA) 
method along the T-L line and Z-F lines are drawn over the spectra in black dotted lines 
and red solid lines, respectively [17]. Comparison of the ARPES spectra to both sets of 
calculations shows that the measured bands do not match either scenario well along lines 
of high symmetry. Interestingly, the shallow bands above -1.0 eV more resemble to those of 
the Z-F line and the deep bands resemble those of the T-L line more closely. 

More useful information about the six valley model is contained in the region of the grey 
rectangular box in Fig. EJ^a), which is seen at length in Fig. [2jb) . Here, three hole pocket- 
like structures are seen as designated with the black dotted lines, I, II, and III, respectively. 
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Fig. 2: (Color on-line) (a) ARPES intensity map obtained at 20 K with ?kj=23.0 eV photon source 
for the valence bands of Te-4% along the FZL mirror plane. The black dotted (red solid) lines are 
the band structure calculated by LDA method along the F-L (Z-F) line [17]. The grey box in the 
top part of the ARPES image is seen at length in (b) . The intensity of the ARPES image is drawn 
in terrain color scale as the scale bar shows, (b) Fine view of the ARPES intensity map for the 
valence band edge. The right two band edges cross the Fermi level while the left one does not. This 
is more clearly seen in the energy distribution curves (EDC) along line I, II, and III. (c) EDC's 
along line I, II, and III. The Fermi edge is seen in EDC II and III. 



In order to check whether those structures are really hole pockets, the energy distribution 
curves along the lines, I, II, and III, are plotted in Fig. [He). Interestingly, band II and III 
cross the Fermi level, which means they form hole pockets, while band I does not as shown 
in Fig. [DJc). These hole pockets, II and III, need comparing with the theoretical prediction 
in detail [17]. According to their LDA result, p-type doped Bi2Te3 has the VBM at (0.546, 
0.383, 0.383) and the second highest band edge at (0.665, 0.586, 0.586) in rhombohedral 
lattice coordinate system. The corresponding points of the VBM and the second VBM are 
located at 0.27 A" 1 and 0.13 A" 1 off the F-Z line towards the Z-U direction with z=0.88Z 
and z=l. 2Z, respectively. The positions of the VBM are depicted as the orange spheres in 
Fig.[Ub). These positions are directly compared to our experimental results. As is described 
above, the two hole pockets, III and II, correspond to the VBM and the second VBM, and 
each fc|| value, ^0.29 A" 1 and ^0.09 A" 1 , respectively, match the theoretical values well. 
Because Bi2Te3 has three- fold rotational symmetry around T-Z axis and two-fold rotational 
symmetry around T-M axis, the system has six VBM's and six second VBM's in the BZ. 
This is direct evidence in support of the six- valley model of p- type doped Bi2Te3. The only 
ambiguity about the six-valley model is which set is the lowest in energy. Unfortunately, 
experimentally answering this question is not trivial because the predicted energy difference 
between them is only 3.8 meV [17], and the time available for measurement of the bulk 
states is very limited. This last constraint is discussed at length below. On theoretical 
aspects also, there is a controversy. Recent LDA calculations have given somewhat different 
results to the question [18,20]. In this respect, our measurements have not been successful 
at proving a definitive model. 

Having established the existence of the multi-valley hole pockets and confirmed the 
effect of the spin-orbit interaction (SOI) on the electronic structure of Bi2Te3, it is desirable 
to understand how the multiple hole pockets induced by the SOI effect enhance the TE 
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Fig. 3: (Color on-line) (a) Time dependence of the ARPES spectra after cleaving. The elapsed time 
is (a) 7 minutes, (b) 30 minutes, and (c) 60 minutes, (d) Constant energy map of the surface state 
at E=Ef. The spectral weight is integrated over an energy range of Ef ± 5 meV. 



properties. Elements of high atomic number like bismuth have large SOI because the spin- 
orbit coupling constant is proportional to the atomic number. As Mahan showed, the TE 
figure of merit, ZT, is a monotonically increasing function of the B parameter: Bj — 
Tefirij{T){kB I e) 2 / '(4/c.t) , %(T) = Nj(2TrmjkBT/h 2 ) 3 / 2 , where e is the carrier charge, \i is 
the mobility, kl is the lattice thermal conductivity, rrij is the density of states mass of the 
j-th band, and Nj is the degeneracy of the band. Increasing Bj is therefore important for 
increasing ZT [23] . Obviously, one of the most efficient ways to increase B is to increase the 
degeneracy, Nj, as in the case of Bi2Te3. The above reasoning and our experimental results 
show that the SOI effect is crucial to understanding the electronic structure of p-type doped 
Bi2Te3 with the multiple hole pockets in the BZ, and is an important factor in explaining 
the high TE figure of merit of Bi2Te3. 

As a phenomenon directly impacting ARPES measurements of Bi2Te3, we also report on 
surface states forming at a cleaved surface of Bi 2 Te 3 crystals. In most cases of our ARPES 
experiments, we observed one or two electron pocket-like structures on the p-type doped 
Bi 2 Te3 crystals, both Bi-1% and Te-4% at any temperature ranging from 10 to 300 K. If 
these states represent the bulk properties of Bi2Te3 their existence definitely contradicts 
our transport measurements. Fortunately, the time-dependent ARPES study gives very 
useful information about the origin of the states. Figure [3fa)-[3jc) show the variation of 
the valence bands of Te-4% in the ARPES spectra with the elapsed time after cleaving the 
samples. When the elapsed time is 7 minutes, it is clearly seen that the valence band edge, 
located off the TZ line, almost reaches to Ep. As time passes, the valence band splits into 
two bands. The valence bands get lowered in energy, and a new band appears crossing 
in about an hour. The Fermi surface of the newly emerging state is circular as shown in 
Fig- Eld). This elapsed time dependence suggests that these bands are not the normal surface 
states induced by translational symmetry breaking owing to cleavage of the bulk crystal. 
Otherwise, they would be seen immediately after cleaving the sample. The endurance test of 
surface contamination also supports this conjecture. If the surface electron pocket is formed 
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by impurities such as excess Bi or Te atoms the impurity concentration estimated from the 
pocket size would be ~1.5%, which is detectable by high resolution x-ray photoemission 
spectroscopy (XPS). In order to check this possibility, we have also performed the shallow 
core level XPS experiments, but did not observe any noticeable change in the Bi 5d and Te 
Ad XPS spectra with increase of the elapsed time. The intensity ratio change of Bi 5d to 
Te Ad peaks was not observed either. This indicates that the states are not relevant to the 
elements in a different chemical environment such as the intercalated excess Bi or Te atoms. 

As to the origin of the surface bands, the surface electronic structure of Bi 2 Te 3 calculated 
by the density functional theory provides relevant information (see Fig. 5 of Ref. [24]). The 
report points out that the surface band structure of Bi2Te3 is different from the bulk one 
due to the weak inter-layer bonding of the layered compound. The main differences are the 
lowering of conduction bands almost to the VBM point, the disappearance of the six-fold 
degeneracy of the conduction band minima (CBM), and non-dispersive behavior of the bands 
along the k z direction. Most features of our ARPES data are consistent with the argument of 
surface effects induced by the weak inter-layer bonding, if we assume the electronic structure 
change in Fig. ETa)-[3](c) is due to slowly relaxing of the quintuple inter-layer distance. This 
assumption is not unreasonable when we consider the plasticity of Bi2Te3. The main reason 
for the CBM lowering its energy in the surface electronic structure is the weakening of 
the covalent ppcr-bonding along the chain direction of Bi2Te3, i. e. Te(I)-Bi-Te(II)-Bi- 
Te(I), induced by an increase of the quintuple interlayer distance. Because the CBM is an 
antibonding state, bond weakening causes the state to be lowered in energy. Inevitably, this 
lifts up the surface Fermi level and induces surface band bending. 

Finally, it is worth mentioning that our finding may require re-interpretation of the 
previously reported photoemission spectra [21]. Our experimental data shows that the true 
bulk electronic structure of Bi2Te3 can be obtained only within ~15 minutes after cleaving 
the samples. If the sample ages more than this time, the structure near the Fermi level 
actually consists of surface states. Further, these states are always seen whenever a good 
cleave was obtained in situ regardless of photon exposure and sample temperatures. The 
energy distribution curve from the ARPES image of the surface states shows very similar 
structure to the conduction bands or the impurity states near the Fermi level in the previous 
data. 

In summary, the electronic structure of p-type doped Bi2Te3 has been studied by high 
resolution ARPES. The valence band edges are located off the T-Z line in the BZ, which 
directly confirms that the spin-orbit interaction is a dominant factor in understanding the 
electronic structure and the corresponding thermoelectric properties of Bi2Te3. Successive 
time dependent ARPES measurements also reveal that the electron-like bands crossing the 
Fermi level near the T point are formed within an hour after cleaving single Bi2Te3 crystals. 
These are interpreted as surface states induced by surface band bending possibly due to 
quintuple inter-layer distance change of Bi2Te3. 
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